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ABSTRACT 

Investigating the cutoff dependence of the Higgs mass triviality bound, the 4 theory is 
formulated on an F± lattice which preserves Lorentz invariance to a higher degree than the 
commonly used hypercubic lattice. I solve this model non-perturbatively by evaluating the 
high temperature expansion through 13th order following the approach of Liischer and Weisz. 
The results are continued across the transition line into the broken phase by integrating the 
perturbative RG equations. In the broken phase, the renormalized coupling never exceeds 
2/3 of the tree level unitarity bound when A/m^ > 2. The results confirm recent Monte Carlo 
data and I obtain as an upper bound for the Higgs mass uir/ f n < 2.46 ± 0.02hte ± 0.08pt 
at A/mR = 2. 
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1 Introduction 



The relevance of the cutoff within the 4 theory arises because of triviality, i.e. the vanishing 
of the self-coupling in the continuum limit. If the model were not trivial, all one would be 
interested in would be the continuum limit, which of course is cutoff independent, and thus 
this investigation into the cutoff dependence of the Higgs mass bound would be irrelevant. 
Because of triviality, the cutoff has to be kept finite in order to retain a meaningful model 
at non-zero coupling. In this way, the Standard Model has to be viewed as an effective 
low energy theory with a built in cutoff, embedded in some yet unknown full theory. For 
self consistency, this cutoff should be at least a few times the Higgs mass, beyond the 
particle content of the effective theory. The closer the Higgs lies to the cutoff scale, the more 
important the effects of the embedding theory (i.e. the cutoff effects) will be. Thus a search 
for a maximal Higgs mass demands a careful investigation of the cutoff dependence on the 
triviality bound. 

To lowest order in the cutoff, the scalar sector is described accurately by the usual <f> 4 
Lagrangian. The effect of the underlying full theory is to introduce corrections in the form of 
(yet unknown) higher dimensional operators. There are two dimension six operators param- 
eterizing the leading order (A -2 ) cutoff effects. The ordinary hypercubic lattice — which is 
commonly used in Monte Carlo investigations of the Higgs bound — breaks Lorentz invari- 
ance already at order A -2 , so the cutoff dependence on the lattice could look very different 
than in the real world. The hypercubic lattice regularization could severely contaminate the 
physics. The _F 4 lattice, in contrast, accurately reproduces all physical 0(hr 2 ) corrections 
and Lorentz breaking effects first occur at order A -4 . Scalar fields can be implemented 
trivially on the F 4 lattice. This idea was first suggested by Neuberger ||T|. 

Despite triviality, the cutoff may turn out to be so low that the physical self-coupling 
could indeed be strong (with respect to a tree level unitarity bound). This calls for non- 
perturbative methods to investigate the triviality bound, even though previous work on a 
hypercubic lattice seems not to leave much room for a strongly coupled scalar sector. Besides 
Monte Carlo methods, where systematic errors arise because of the finite system size, the 
strong coupling (or high-temperature) expansion has been worked out. Here results are 
obtained directly in the infinite volume limit, without any statistical errors. The method, 
however, only works in the symmetric phase and relies on the perturbative renormalization 
group to be continued into the Goldstone phase. It is also quite cumbersome to apply the 
high-temperature expansion to lattice actions that include additional higher dimensional 
operators which modify the naive nearest neighbor interaction. 

Monte Carlo simulations have been previously performed on an F4 lattice, but a strong 
coupling expansion of the model on this particular lattice has been lacking so far. I will 
essentially follow the approach of Liischer and Weisz, presented in a series of papers ]2|,|3],f|,||] 
and extend the application to the F 4 lattice. For completeness, I will repeat some of their 
techniques and emphasize the adaptions necessary for the F4 lattice. 

I work within the Dashen-Neuberger approximation ||, i.e. Gauge and Yukawa couplings 
are neglected. They are weak and may be included later via perturbation theory. This 
approximation holds (i.e. captures the relevant physics) if the top quark is well below 
0.5 TeV (which seems to be favored by experiment). In this limit the Higgs action turns into 
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the action for an 0(4)-symmetric 4 model 

Sh = fd^^d^-d^+^-v 2 ) 2 ] . (1) 

The two physical scales in the broken phase are set by the Higgs mass and the pion decay 
constant f n = 246 GeV. Based on f n , a physical self-coupling qr is defined as 

if) 2 - 

Triviality of (J) 4 in four dimensions is not rigorously proven but rests on perturbative renormal- 
ization group (RG) results and Monte Carlo data f7|,|8|,p,iT|li~l|. Perturbative renormalization 
group leads to the relation 

mR /A = C 1 (P 1 g R y^exp(--^—){l + 0(g R )} as g R - , (3) 

\ Pi9r ) 

with C\ being a A-dependent integration constant that will be evaluated in the high-tempera- 
ture expansion. Removal of the cutoff A will lead to a trivial theory. Lowering the cutoff, 
however, g R is increasing accordingly. But so are cutoff effects and by putting a limit on the 
acceptable cutoff dependence in physical observables (for instance pion-pion scattering at 90 
degrees) one gets an upper bound on g R and hence on tur/ f n defined in a lattice independent 
way, valid beyond the realms of lattice gauge theory. 

The smaller the coefficient C\ in (|3|) (which depends on the higher dimensional operators 
in an improved action) the less the impact a change in the cutoff will have on mR/f n . This 



in turn allows a higher mass bound. I refer to ||r2],p^,|14|| for a detailed analysis. 

This paper is organized as follows. In the next section I will briefly review some pertur- 
bative results for the </> 4 model on an F4 lattice. I will give scaling laws in the symmetric and 
in the broken phase. (Expressions up to three loops for the RG functions are listed in the 
appendix.) These results were worked out by Luscher and Weisz p],|3|,f|,|5| and by Bhanot et 
al. fll2"l) |i~3| . In the following section, I will give a detailed account of the high-temperature 
series in form of a linked cluster expansion. The method follows work by Luscher and Weisz 
(henceforth referred to as LW) [§], but is adapted to the F 4 lattice. The analysis of the 
expansion up to 13th order is presented in the last section (preliminary results were given 
in [|i~5l, |i~6|1 ) . Results in the broken phase are obtained by propagating the high-temperature 
expansion results using the perturbative RG equations given in section 2. I will show that 
the region where the cluster expansion converges and the perturbative region do overlap on 
the F 4 lattice. An estimate for the Higgs triviality bound is given and compared with Monte 
Carlo data. Ultimately, this comparison serves as an indicator of the size of systematic er- 
rors. Furthermore, the F4 result shows that Lorentz invariance breaking effects are less than 
10%, indicating the robustness of the triviality bound. 



2 Definitions, scaling laws and matching 



Many of the perturbative results required for the high-temperature expansion analysis listed 
throughout this section are taken from || and |I2| and are reproduced here for completeness 
only. 
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The naive euclidean lattice action for the 4 theory is 

S = -2k £ ct>{x)cj)(x') +Y / uU(x)) , (4) 

<xx'> x ^ ' 

with a local 0(N) invariant potential 

u{<j>) = <P 2 + A(0 2 - l) 2 . (5) 

The summation in fl|) is over all nearest neighbor pairs < xx' > and is a real iV-component 
field 4> a (x), a = 1,...,N, located on an lattice. The parameter range is restricted to 
k > and A > 0. The model is known to exist in two phases, separated by a second order 
critical line k c (X). 

The lattice in d > 2 dimensions is defined as the set of points 

{x\x = X)u a; M e (Ai) ' x m e S» x m = even}. Here e^) is the euclidean unit vector in the 
//-direction. F d can be viewed as a hypercubic lattice Z d with its odd sites removed and for 
d = 3 it corresponds to an fee lattice. The euclidian distance between nearest neighbors is 
Vz and there are 2d(d — 1) nearest neighbors per site. The discrete rotation symmetry group 
is particularly large in four dimensions and in this case the Lorentz invariance breaking term 
in the kinetic energy part of the free propagator vanishes. Invariance breaking terms first 
occur at order A -4 . 

The generating functional W for the connected correlation functions of <p and O in the 
presence of the external sources J(x) and K(x) is 



x) , (6) 



exp(V[J,K]) = J n^(x)exp|-5[0] + £ J'W(s) 

where O is a composite operator defined as 0{x) = 2J2<xx>>4 > ( x )4 ) ( x ')- The connected 
correlation functions are then obtained by taking the derivatives of W[J, K] with respect to 
the sources 



W^ l \x u ...,x 



ni Vli ■ ■ ■ j yi)a\...a r , 
6 n+l 



-W[J,K] 



(7) 

J=K=0 



5J ai (xi) . ..5J an (x n )5K(y 1 ) . ..SK(y t 

The vertex functions 1^"'^ are generated by the one particle irreducible functional T[M, K\ 
defined as the Legendre transform V = W — J2 X J{x)M(x). Then 



; 9l> • • • j Ql)ai...a n — 
5 n+l 

T[M, K] 



M=K=0 



6M( Pl ) ai . . . SM( Pn ) an 5K( qi ) . . . 5K( qi ) 

and M is the local magnetization, M = dW/dJ. Because of O(N) invariance, the coefficients 
p(n,o van i s h f or oc id n. 

The renormalization conditions imposed on the vertex functions in the symmetric phase 
at zero external momentum define the wavefunction renormalization constant Zr and the 
renormalized mass tjir, 

r {2 ' 0) (p, -PU = 5 a ^{p 2 + m R } + 0(p 4 ) , (p - 0) . (9) 
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The renormalized coupling is fixed via 



r( 4 <°)(0, 0, 0, 0U, S = \c 4 (a, (3, 7 , 8)ik > ( 10 ) 



with C4 being the totally symmetric O(N) invariant tensor 64(01, (3, 7, 5) = 5ui3$~i& + Sa-ydps + 
8as8/3-y The field renormalization constant Z® associated with the composite operator O is 
given by 

(0,0; 0)^ = 5^^^. (11) 
The renormalized vertex functions are defined by 

r£'° = for odd n and for n = 0, / < 1 , 

r W) = ^/2 z o^ r (n, _ 5nofer (o,2) (0!0) | otherwise) (12 ) 

and satisfy the following normalization conditions 

rg ,0) (p,-p)^ = M? 2 + ™*> + S(P 4 ) > (p->0), (13) 
rg' 0) (0,0,0,0) a/37<5 = ^C 4 (a,/3, 7 ,%Ji, (14) 

rf 2) (o,o) = 0, (15) 
rg-^cojo)^ = 5 a/3 . (16) 

From the lattice action, the vertex functions will be expressed in terms of k and A. Taking 
the derivative of the renormalized vertex functions with respect to k at a fixed value of A, 
one derives the Callan-Symanzik equation 



om R dg R 



(n,I) 



,2 



|rt' +1) n -Wi 2 rg' s) (o,o,o)} . (17) 



The coefficients are defined by 



(3(m R ,g R ) = ™«|f/^f , (18) 



1 d\nZ R /dm R 
l(m R ,g R ) = - mR -^- — , (19) 



dlnZp I dm R , , 

6(m*,gn) = m R -^ / ' , (20) 

e(m R ,g R ) = L R Z°?j£\ = 2( 7 - 1) . (21) 
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Equivalently, one may write the above equations in the following way 



miiP^ = P, (22) 



"dm 



R 



d\nZ R , . 

m R~K = 2 7 , 23 
om R 

mR ~PtZT^ = ^ (24) 



dm R 
dn 



dm 



R 



m\eZ° R , (25) 



where all derivatives are taken at fixed A and m R is the independent variable. At the 
transition, as m R — > 0, the RG functions /3,j,8 and e have finite limits. The universal 
coefficients P u , r y u ,S u and e u have been previously calculated up to three loops. They are 
given in the appendix together with the non-universal scaling violating terms up to one 
loop. 

Equation (^) implies that, for a positive /3-function, as m R — > the coupling also 
tends to zero as given in eq. (|3[). Similarly, scaling laws for Z R , Z R and k are derived from 
eqs. ([£f)-(]2"5j). Specifically, one has 

Z R = C 2 {l + 0(g R )}, (26) 
Zg = C 3 (g R ) s ^{l + 0(g R )}, (27) 
k c -k = C 3 m 2 R (g R y^{l + 0(g R )}. (28) 

Switching from m R to a new independent variable r = 1—k/k, c , the equations are transformed 
to 



m R ~ 



C,T l ' 2 \\nT\ & ^ , (29) 



g R r ~° J-llnrl" 1 , (30) 
Pi 



7 r^O 
'R ~ 



C 2 , (31) 
Z% r ~° C 5 \lnr\- 5 ^ . (32) 

At this point, note that the susceptibility (as defined in ( |4*6"D below) is related to the renor- 
malized quantities via \2 = Z R m R 2 , resulting in the following scaling behavior 



X2 



T ~° C 6 r- 1 |lnr|- <5l//31 . (33) 



At small A however, the divergence of the susceptibility is dominated by the r _1 term, unless 
one is very close to the critical line. This can be seen as follows. Keeping terms proportional 
to the initial value go = 6X/(d — 1) 2 k, 2 when integrating the RG equations, one obtains 

9R ~° Soll + ^llnrir , (34) 
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and thus the improved scaling law for the susceptibility 



X2 



C-jt- 1 l-/(A)|lnr| 



-61/ Pi 



Here I have defined 



/(A) 



d2N + 8 x + o(x>) 



(35) 



(36) 



An 2 N + 2 

where A = (N + 2)A + 0(X 2 ) (see ( f79|) for a proper definition) and I used the fact that k is 
close to k c = l/2d(d — 1) + 0{\). All the other scaling laws are modified accordingly. The 
improvement indicates that the logarithmic divergence is suppressed for small A. 

In the broken symmetry phase, the iV-th component of <fi is given a vacuum expectation 
value v while the first N — 1 components remain massless. The generating functionals are 
defined in the limit of a vanishing external magnetic field h. The vertex functions r*™'^ 
are again given by (|8|). However, in the broken phase they are only 0(N — 1) invariant. 
The perturbative expansion of the vertex functions in the broken phase is defined around 
4>a = 8aNSmin, where s m i n is one of the two degenerate absolute minima of the lattice action. 

The renormalization conditions in the broken phase are imposed (following |12|) as 



rg ,0) (p, -p) aP 

Re J r ( 2 '°)( P , - P ) NN 



5 apP 2 + 0(p 4 ), (p-0) 



Re|r (1 ' 1} (p; -p) 



N 



p=(im R , 0,0,0) 
p=(im R , 0,0,0) 



ZrZr 



(37) 
(38) 

(39) 



The pion decay constant is just the vacuum expectation value of the renormalized cr-field 
f n = (Zr)~^v. Now the renormalized self-coupling is defined as in eq. (0) and renormalized 
vertex functions are introduced via 



R 

,(n,0 
R 



for n = 0, I < 1 , 

Z ^ Z Ol[ r (n,l)_ Sno5[2rm{p _ p) 



(40) 



p=(im H ,0,0,0) 

The Callan-Symanzik equation in the broken phase reads 



otherwise . 







m R 







dm R dg R 

= §m R ri n+1, ° 



15\Y 



[n,l) 
R 



p n+1 =0,a n+1 =N 



1 ^ "nV " 
+ ^R T R 



(n,l+l) 



1l+l=0 



(41) 



The RG functions are defined as in eqs. (0) — (|21|). However, the relation between e and 7 
is not valid in the broken phase, and — unfortunately — in the case of iV = 4 there exists 
no simple relation between e and the other RG functions. The function $ is defined as 



^{m R ,g R ) 



-1/2 dv_ / dm R 
1 0k / 8k 



(42) 
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and it is related to the other functions through 




# = 1 + 7-^h/—. (43) 



The universal coefficients P u , r y u ,S u and e u plus the scaling violating terms are listed in the 
appendix. 

The scaling laws for tur, Zr, and Z R — given in (|3]), ( p6| ) and (|27| ) — are exactly the 
same, except with integration constants C- instead of Cj. The only change is that in the 
scaling law (ESI) for k appears an additional factor of —1/2 



«-«c = ^C' 3 m 2 R (g R ) s ^{l + 0(g R )}. (44) 

As in the symmetric phase, the scaling law for m R implies that the renormalized coupling 
tends to zero when approaching the critical line from above. 

To connect the theory in the symmetric with the theory in the broken phase, one has 
to establish a relationship between the integration constants Ci and C[ on both sides of the 
transition. This is accomplished by reconstructing the massive renormalized vertex functions 
in both phases from the critical renormalized vertex functions by mass perturbation theory. 
The answer can be obtained from a one-loop calculation. One finds 

= Cl (A)exp{ 2jV ^ + -f- }, ,45) 

and C[ = d for i = 2, 3. 

Finally, the tree level unitarity bound in the symmetric phase from the S'-wave phase 
shift for elastic scattering amounts to gn < 29 for N = 4. An almost equal bound holds in 
the broken phase f|. 



3 The linked cluster expansion 

As a powerful non-perturbative semi-analytical tool, the high-temperature expansion serves 
as an alternative to Monte Carlo methods. It has been widely used and for details I refer 
to P, [17|. The high-temperature expansion can be systematically organized in terms of 
graphs consisting of vertices and connecting lines. This method is referred to as linked 
cluster expansion. The number of lines in a graph corresponds to the power of the hopping 
parameter k in the high-temperature expansion. I compute susceptibilities up to 13th order 
in n with a total of about 400 000 contributing graphs. The number of contributing graphs 
depends on the lattice structure as not all graphs can be embedded on all lattices. The whole 
procedure is computerized as outlined in || and is easily adaptable for different potentials 
and lattice structures. 

It is sufficient to define the following susceptibilities in terms of the connected correlation 
functions and 

X2 = ^£(W<r(0)> com \ (46) 

x,a 
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to = ^E^(^(^)^(0)) c ° nn , (47) 

x,a 

X4 = ^Y.{n^ a {y)<t>\z)<t> b m^\ (48) 



a,b 



where the spatial summation produces zero momentum operators. From the expansion of 
the susceptibilities an expansion for the quantities of interest, particularly the mass tjir, the 
coupling qr, and the wavefunction renormalization constants Zr and Z® can be calculated 
in principal anywhere in the symmetric phase via 




m R = , (49) 

g R = -(2d) 2 |i, (50) 
Z R = 2d^, (51) 



These relations can easily be derived from the renormalization conditions imposed on the 
vertex functions in the symmetric phase. To extend the solution into the broken phase one 
has to resort to other means, e.g. perturbative RG techniques. 

Care has to be taken in defining a proper Fourier transform on the underlying lattice 
structure. The particular form of the Brillouin zone on the _F 4 lattice results in additional 
factors of 2 in the above relations as compared to the hypercubic case. 

I will now introduce the necessary ingredients for the cluster expansion, starting with 
some basic definitions. A graph G is defined by a set (V, £,£), where V represents the 
vertices, C the internal lines connecting the vertices and S the number of external lines 
attached to each vertex. Let V be the total number of vertices, L the number of internal 
lines or order of the graph, E the total number of external lines. Let N(v) = I(v) + E(v) be 
the total number of lines, internal plus external, attached to a vertex v. Two graphs G and 
G' are called topologically equivalent, if G' can be transformed into G by simply reordering 
its vertices and [G] is the corresponding equivalence class of graphs. In the linked cluster 
expansion, the particular order of the vertices is irrelevant and therefore the expansion is 
in terms of equivalent classes rather than individual graphs. Care has to be taken to avoid 
overcounting. 

The distance between two vertices d(vi, Vj) is the length of the shortest path between 
Vi and Vj. A loop is a closed path. Generally, a graph may contain vertices that are not 
connected by a path to some (or all) other vertices. A graph G is said to be connected if 
there is at least one path between any two vertices. A connected graph G is said to be 
non-separable if for all vertices v the graph G remains connected after removal of v and all 
lines emanating from it. A graph is called simple if any two vertices are connected by at most 
one line and there are no self-loops (i.e. no lines starting and ending at the same vertex). 

The linked cluster expansion is based solely on graphs without self-loops, and therefore 
contributing graphs may be represented by a V x V symmetric matrix. The entries A(i,j) 
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of this so-called incidence matrix give the number of internal lines from to Vj for i j, 
and A(i,i) = E(vi) otherwise. The incidence matrices of topologically equivalent graphs 
are related through a permutation of the vertex labelling (1, . . . , V). To eliminate multi- 
ple occurrences of equivalent graphs in the expansion, one has to define a certain canonical 
ordering of the vertices which must depend only on the topology of the graph under consid- 
eration. As a result, every equivalence class of graphs is represented by a single incidence 
matrix A(i,j). 

With each graph in the expansion one associates certain characteristic numbers which 
may depend on three sources: the topology of the graph, the underlying lattice structure or 
the local potential in the action. 

The potential u{4>) only enters through the following weight factors. For a graph G with 
a set of vertices V = (v\, . . . , ty) it is defined as 

° , , -w-r o conn 

W (G) = J] m m , (53) 

o conn 

where the m 2k are the so-called connected moments. They can be calculated from the 
potential as follows. 

The Taylor expansion of the free energy W[J,K; k] at zero coupling where all spins are 
uncorrelated and essentially randomly distributed can be written in terms of the so-called 

o oi.-.aj. 

cumulants M k (xi) defined via 



d k W 



dJ ai (xi) . ..dJ a *(x k ) 



= 6(xi,...,x k ) M,! h (xi) . (54) 

K=0 



For zero external field, only the even correlation functions W^ 2k ' ^ survive in a symmetric 
potential and the cumulants take a simple form 



M 2k (xi) 



= (<p(^---<P fc (^)> 

J=K=0 



Ik=0 

<s(zi, . . . ,x 2fc )(<r (*i) . . .r 2k (xi)r nn , (55) 



where the connected correlation functions are evaluated at one point in space-time and may 
be denoted by (0 ai . . . a 2fe^onn^ p rom O(iV) invariance it follows that 



o conn 



(<r...0 a -) c r n = C 2k (a 1 ,...,a 2k )m 2k 5( Xl , . . . , x 2k ) , (56) 

o conn 

where m 2k are the connected moments appearing in (pjj). The C 2k (ai, . . . ,a 2k ) are O(N) 
invariant tensors. The connected moments can easily be calculated from the (not necessarily 
connected) moments, which are defined in terms of the correlation functions 

(0 ai . . . r 2k )i =C 2k ( ai ,..., a 2k ) m 2k . (57) 
Using the proper normalization condition for the O(N) tensors one derives 

o r(±N) J N -i+2k , _ v 

m2k ~ ¥rWTk) j n _, ' (58) 
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with 

J t = J™ dr r l exp (-u(r)) . (59) 

In the limit A — > oo, the moments can be written as a ratio of polynomials in N. One has 
exp(— u((f>)) oc 5(0 2 — 1) which results in J; = const, and thus 

o 1 , , 

m2k ~ (N + 2k - 2)(N + 2k - 4) . . . (N + 2)N ' ^ ' 

From the moments the connected moments are generated via the relation 

~ o conn Z k , f ~ Z k ) 

L m 2fc Tf = ln l + ^m 2fe - . (61) 
fc=i ' I fc=i ' J 

This concludes the calculation of the weight factors. 

Next, with every graph G, one defines an internal symmetry number S(G) as 

S(G)=S v (A)UMhjy- > (62) 

where S^A) is the number of permutations that leave the incidence matrix invariant and the 
product represents the internal line multiplicity factor. A modified symmetry number Se{G) 
is defined which additionally takes the external line multiplicity into account by setting i < j 
in the above expression. 

The underlying lattice structure enters through the lattice embedding number 1(G) which 
counts the number of embeddings of G on the lattice under these conditions: the first vertex 
in G is fixed at the origin of the lattice and all other vertices are then distributed in such a way 
that any two vertices with distance 1 are placed on nearest neighbor lattice sites. Different 
vertices may be placed on the same site and different lines may occupy the same lattice link. 
Not all graphs are possible on a certain lattice structure. For example, graphs containing 
odd loops cannot be embedded on a hypercubic lattice. Note that for the susceptibility fi2, 
one has to introduce a modified lattice embedding number h^G) which takes into account 
the additional factors of x 2 . 

Finally, one defines a symmetry number associated with the internal 0(N) symmetry. 
For a graph G with a set of internal lines C = (h, . . . one writes down an O(N) index 
ai G {1,2,..., iV} for every line /. Let {ii,i2, ■ ■ ■ , h(vi)} be the set of indices labelling the 
internal lines which terminate in vertex V{. Then the O(N) factor for G is given by 

N V 

C{G) = II ^Njvj) ( M, • • • , 1 , Qii , Qj 2 , • • • , Qj JK) ) ■ (63) 

a 1: ...,a L =l i=l ^~ ^~ ' 

hi(Vi) times 

C(G) is a polynomial in iV with maximum degree \L. 

Having defined all the basic ingredients, one may now write down the unrenormalized 
expansion for the susceptibilities based on all equivalence classes of connected graphs without 
self-loops and with N(v) even for all v E V and E = 2k. These graphs are called Q 2 k graphs 
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and the expansion reads 



X2 = J2(^) L HG)C(G)W(G)-^- (64) 
Geg 2 z>e{G) 

/x 2 = E (2*;) L / 2 (G)C7(G) # (G)-^- , (65) 

X4 = E (2*0 L /(G)C(G) # (G)^y • (66) 



Geg 4 



It can be recast in the form 



X2 = 2^X2 , — F r- , (67) 

L=0 ^ ! 

= X>?^ , (68) 

L=0 ^ ! 
L=0 ^ ! 

where the coefficients xi^ > an d xi^ are °f the form 

oo 

p^in^rv, (70) 



with J2kLi k n k = L + i-E. The p[iV] are polynomials in N of order not higher than \L. 

In order to reduce the total number of graphs, two levels of renormalization are carried 
out. A reduction in the number of contributing diagrams in the expansion is paid for with 
an increase in algebraic complexity. 

The first step is to eliminate one-particle reducible graphs. A graph G G Q 2 k is called 
one-particle irreducible, or 1PI, if it cannot be broken into two disconnected parts by cutting 
a single internal line. The set of 1PI graphs is denoted by Q\^ 1 ' . The expansion over Q 2 k can 
be reconstructed via 

X2 = xril-ex 1 / 1 )' 1 , (71) 



A*2 = 



1PI i ^(^,IPI\2 



/4 + <x 



2 



exV 1 )- 2 , (72) 



X4 = xr(l-%^)- 4 , (73) 

where x\ PI ■> tA PI an d x\ PI stand for the 1PI part of the corresponding reducible quantities 
and e = qln. Here q is the coordination number of the lattice. 

In a second step, one throws out one-vertex reducible graphs. A graph G G Q^ 1 is called 
one-vertex irreducible, or 1VI, if — after removal of any vertex and all the lines emanating 
from it — each remaining connected piece of G contains at least one external line. The set 
of 1VI graphs is denoted by S 2 k- Given the set <S 2 fc, the set can be reconstructed by 
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attaching <S 2 fe graphs repeatedly to all G G S2k in a certain manner which I will not describe 

o conn 

here. Furthermore, the connected moments m 2k for the S2k graphs have to be redefined as 



m2fc = Jot 1 m E (^) L HG)C(G) w(G)-±- (74) 

The rri2k are called renormalized moments and are evaluated over the set V2k, which is given 
by all graphs G G Q2k with the additional constraint that all 2k external lines are attached to 

o conn 

a single vertex. To lowest order one has rri2k =m 2k +0(k ). The structure for the moments 
is the same as for the susceptibilities above. One finally has 



m 

Xl pi = E (2k) l I(G)C(G)W(G)— — , (75) 

H\ PI = E (2k) l I 2 (G)C(G)W(G)-^- , (76) 

FA 

XT = E {2k) l I{G)C{G)W{G) — — . (77) 

Ges 4 dE{<~r) 



I will not go into details of the coding; the interested reader is referred to and [TjJ. I 
will only give a brief outline and stress the differences to the hypercubic case. 

To construct the set of S2k graphs required for the IPI susceptibilities, I start out with a 
set of 'simple' graphs which will be expanded into the full set by a couple of basic operations. 
Let C\ be the set of simple and connected graphs without external lines which are IPI and 
non-separable. As external lines are absent, non-separability is equivalent to one- vertex 
irreducibility. I denote the subset of C\ graphs with L internal lines by Ci(L). All the 
C\(L) graphs can now be constructed by repeatedly applying two simple operations to the 
triangular graph, Ci(3), represented by 



( 


1 


i 


1 





i 


V i 


i 






(The zeroth order graph has to be added by hand.) All higher order graphs can then be 
generated 

• by adding a single line between any two vertices with A(i,j) = 0, and 

• by replacing a line by a bridge: a new graph is generated by removal of a line, say from 
Vi to Vj and addition of a new vertex w which is connected to both, Vi and Vj, by a 
single line each. The number of lines and the number of vertices is thus increased by 
one. 

It is obvious that the resulting graphs are still simple and connected. However, the same 
graph may be constructed more than once. On the other hand no graph in C\ is left out, 
as any graph may always be reduced to the basic third order graph by simply applying the 
inverse operations. On a hypercubic lattice, graphs with odd loops are ruled out and one 
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would therefore start with the fourth order square graph. In this case it is convenient to 
modify the rules of construction such that graphs with zero lattice embedding number are 
not constructed at all. The number of C\(L) graphs up to 13th order is given in table 1 for 
the general (i.e. F4) case and for the hypercubic lattice. 

To construct the set S 2 k from the set Ci, two additional conditions have to be satisfied. 
The total number of lines attached to each vertex N(v) has to be even, and the number 
of external lines has to be 2k. For the graph to be 1VI, at least two vertices must have 
external lines. Besides adding external and multiple internal lines to fulfill these conditions, 
the resulting graphs need to be coupled to generate the complete set of S 2 k graphs which are 
not necessarily non-separable but only 1VI. Details are given in ||15|| . The number of S 2 k for 
k = 1, 2 is also shown in table 1. 

The total number of graphs generated up to 13th order is about 400,000. With 32 bit 
precision, this would require \ ■ 13 ■ 14 ■ 4 ■ 10 5 ■ 4 bytes = 145.6 Mb of storage on the computer, 
However, the incidence matrices can be stored in compressed form. Each matrix element 
requires a maximum of 4 bits at 13th order. Thus 12 4-byte integers are sufficient to represent 
a single matrix. A subset of graphs is simple and can be represented by just 3 integers. The 
total storage is therefore less than 19.2Mb. 

The internal symmetry numbers and the O(N) symmetry numbers are independent of the 
lattice structure and are calculated just as in the hypercubic case. For a detailed description 
of the procedure see 0. The lattice embedding numbers are more involved and depend on 
the choice of lattice. They are calculated for a subset of so-called reduced graphs in order 
to minimize the calculational effort. There is a total of 6180 of these graphs up to 13th 
order, as compared with a total of 396,140 S 2 k graphs (k = 1,2). The evaluation of the 
embedding number involves the number of random walks Mi[x\ from the lattice origin to x 
in I steps. Instead of using a combinatorial approach to calculate these numbers, they are 
computed by actually carrying out the random walks on the lattice and creating a look-up 
table. This method can be adapted for different lattice structures without much effort. For 
the modified lattice embedding numbers each contribution is multiplied by an additional 
factor of (x(i) — x(j)) 2 , where u$ and Vj are the vertices with external lines. Factors of \[2 
enter on the F4 lattice due to the length of basis vectors. 

The renormalized moments are evaluated as in ([74]) over the set V 2 k which unfortunately 
is one-vertex reducible and therefore quite large. By introducing yet another set of 1VI 
graphs and new moments associated with these, it is possible to compute the renormalized 
moments order by order in terms of the new moments in a recursive manner. Consider a 
graph G E Q 2 with E(w) = 2 and I(w) = 2k for some vertex w. Now construct a new graph 
G' by removing the vertex w plus both its external lines. Note that G' is not necessarily 
connected nor 1PI or 1VI. One then defines 



Q 2k = {G E G 2 \G' connected and 1VI } . (78) 

G' does not need to be 1PI. By construction, all graphs in Q 2 k{L) will satisfy L > 2k. 
Observe that for G' to be one-vertex irreducible, all spikes must end in an external line. 
All external lines must have been connected to the vertex w before its removal. Therefore, 
G will consist of loops only. This in turn means that all Q 2 k graphs will in fact be 1PI. 
For a hypercubic lattice, all loops have an even number of lines and since N(v) is even, 
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Table 1: Number of graphs contributing to C±, S2, S4 and 0,2k- The second row gives the 
number of graphs that can be embedded on a hypercubic lattice. The numbers for Q,2k are 
the sum of all contributing graphs, k — 1, . . . , 00, at the given order. 
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Figure 1: The contributions to 5 2 and 5 4 up to 5th order. 
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too, the graphs G, as well as G' must be of even order. This condition does not necessarily 
hold on other lattices, though, and in particular it is not true on an F 4 lattice. Finally, 
the new moments are evaluated in terms of the Q.2k(L) graphs and these are related to the 
renormalized moments in a recursive manner. Again, the full procedure is described in detail 
in@. 

With the renormalized moments, I have all the necessary components to evaluate the one 
particle irreducible parts of the susceptibilities as given in eqs. (|75|) — Q77p. After a simple 
re-expansion of the series, I obtain the final answer in the form of eqs. (|67j) — flBTf). Limited 
by computational resources, I carry out the expansion up to 13th order. On a hypercubic 
lattice, where the number of graphs is significantly fewer, Luscher and Weisz were able to 
go to 14th order. 

The total number of contributing terms in the expansion for L < 13 on the F4 lattice is 
387 for X2, 357 for /i 2 , and 562 for Xi- The coefficients are given in tables 2 — 4 for the non- 
linear limit. For general A, the coefficients are defined in terms of polynomials p[N] multiplied 
by powers of the connected moments, see eq. (|70|) . These numbers can be obtained from the 
author by electronic mail. 

As a check on the computer program, I reproduce the high-temperature expansion up to 
14th order on a hypercubic lattice. I am in full agreement with Luscher and Weisz 0. I 
also check the expansion on an F3 (or fee) lattice, where most previous data is for the Ising 



model. For X2, I find agreement up to L = 12 with Gaunt, Sykes et. al [|TJ]. For both, X2 
and H2, I agree with the data in the Baker-Kincaid tables [19j for the N = 1 case up to 10th 
order. My results are also in agreement with the data by Moore for the Ising model on an 



-F4 (or hfee) lattice up to 10th order [20]. For fi2, however, I disagree at 10th order in the 
14th digit. Since my d = 3 results for [12 agree with Baker-Kincaid at that order, and the 
d = 4 case is just a trivial extension, I am confident that my results are correct. 



4 The triviality bound 

From the high-temperature expansion the relevant physical quantities are calculated follow- 
ing Luscher and Weisz [§, Some modifications are necessary for the F4 lattice. In the 
following analysis I will only consider the physical four component case. The evaluation 
of the series is improved by taking into account the scaling laws quoted in section 2. The 
results from the high-temperature expansion in the symmetric phase are being used as initial 
data to perform an integration of the RG equations towards the critical line. By matching 
the scaling behavior of the renormalized coupling, the integration is then continued into the 
broken phase up to m# = 0.5. Finally, the result for mu/f^ is compared with Monte Carlo 
data. 

4.1 High-temperature series analysis in the symmetric phase 

Approaching the second order critical line k c (X), the correlation length diverges - - i.e. 
A/mR — > 00. In this limit, the theory is expected to behave like an effective, low-energy 
continuum theory. As one wants to obtain an upper bound on A/m# at fixed renormalized 
coupling gji (i.e. fixed physics) one has to minimize the function m^(A) along the lines of 
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Table 2: Coefficients X2 for d = 4, A = oo and N = 4 on an F 4 lattice. 
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Table 3: Coefficients /4 for d = 4, A = 00 and N = 4 on an F4 lattice. 
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constant g R (the RG trajectories). Thus, the first step is to compute the functions m R (n, A) 
and g R {K, A) from the high-temperature series. 

The procedure is as follows. First I will determine the critical line k c {X). To identify 
possible additional singularities I apply the Pade approximant method. (Remember that on 
the F4 lattice there is no simple anti- ferromagnetic singularity at k = —k c .) Closer to the 
critical line, convergence of the series expansion becomes rather slow. However, LW found 
that the rate of convergence is "acceptable" as long as the correlation length does not exceed 
the average size (as measured in lattice spacings) of the graphs contributing to the cluster 
expansion. I choose as limit m R = 0.20. In terms of k, this condition is always satisfied 
in the region k < 0.98k c . To obtain results in the region above 0.98k c , I integrate the RG 
equations (again at fixed A) . The starting point for the integration is taken at the boundary 
line k = 0.98k c . As g R is driven towards zero at the transition, the perturbative expansions 
of the RG functions remain valid — provided the coupling is weak enough at the boundary 
line. A crude estimate of the applicability of renormalized perturbation theory is given by 
the tree-level unitarity bound on the renormalized coupling. 

At this point let me recall that the results from the high-temperature expansion are for 
the susceptibilities X2, ^2 and \a as m e Q s - (EH) — (|48|) . From these, m R , g R , Z R , and Zf 



J R 

can be computed via eqs. (HDI) — fl52"|). For the analysis, I switch to a new bare coupling A 



o conn 

which is defined in terms of the connected moments m 2k (see eq. (|56|) ) as 

— N o conn ; o conn „ 

A = -(1 + Y )m 4 /(m 2 ) 2 • (79) 

It runs monotonically from to 1. See table 5 for some values. 

The convergence properties of the high-temperature expansion have been found to im- 
prove by reexpanding the series in powers of a so-called character variable v B. It can be 
expressed as 



N+l 



^ /' V ( N+l+l \ j2 
Z=l,3,5,... t- 1 V 2 / 



1 ,r(i±i) 2 (80) 

1=0,2,4,... 1 1 2 / 

In the Ising limit this reduces to the well known transformation v = tanh2/t. k c (A) is given 
by the pole in the series expansion of the susceptibility X2 and is extracted using the ratio 
method. Because all coefficients in the expansion X2 = X^o X^l v% are positive, the pole 
closest to the origin must be on the positive real axis. This implies that the singularity is 
given by v c = linij^ooTj with r { = x£ /xiP- For i = 13, the ratios have already converged 
rather well. 

Based on the knowledge about the scaling behavior of the physical singularity in X2, the 
extrapolation i — > 00 can be improved. This behavior is described by eq. (|35J) and can be 
simulated by an auxiliary function of the form 

, / \ -61/ Pi 

h(z) = (1 - z) (l - /(A) ln(l - zyj (81) 

00 

i=0 
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Table 4: Coefficients x<L for d = 4, A = oo and N = 4 on an F 4 lattice. 
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Table 5: Values for k c and A at given A. /cf is from the Pade approximant, k c is from the 
extrapolation of the ratios. 
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where /(A) is given in eq. (P^). Taking the ratios r[ = h^~^ /h®, one would expect similar 
convergence properties for and r[. These two sequences are then combined to form a ratio 
R% — r i/ r i, which somewhat stabilizes the series (especially closer to the non-linear limit). 
The value for v c is obtained by fitting the monotonic fall-off in Ri to a 1/i 2 form. The error 
Av c is taken as i\Ri3 — v c \. 

To study the distribution of singularities of \2 i n the complex plane I use the Pade 
approximant method. Suppose the singularity of f(z) can be written in the form f(z) = 
A (z — z c ) a {1 + 0(z — z c )} for z ~ z c . Then the logarithmic derivative of f(z) exhibits a 
simple pole at z = z c which is extracted using the Pade approximant. The exponent a is 
given by the residue at the pole. 

Results for k c from both methods at selected values of A are displayed in table 5. Notice 
that the values for k c from the Pade approximation are consistently below the result from 
the extrapolation of the ratios Ri. The difference is largest for intermediate A where the first 
order scaling law improvement may not be adequate. On the other hand, the Pade method 
completely ignores the additional information about the nature of the physical poles gained 
from the scaling laws. Furthermore, the error estimation for k£ is somewhat unclear and it 
is therefore difficult to draw definite conclusions. I will thus consider the results from the 
ratio method superior and use them throughout. 

At A = oo, the critical point is k c = 0.09155(5) which agrees with previous Monte Carlo 
data, k^ ic = 0.0917(2) |0|]. Except for the pole at k^, the Pade approximants only give 
spurious zeros or poles far away from the origin. Almost all of the roots in the polynomial 
d(z) are on or very close to the real axis. The value for a from the Pade approximant is 
close to -1, as expected. 

The calculation of the renormalized quantities rriR, g R , Zr, and Z^ is carried out along 
the boundary line k = 0.98k c . From the scaling laws eqs. (|29| ) — (p2|), the expansions can be 
written as 

\ 12 v~ 1/2 m R {v) , (82) 

(83) 

) , (84) 

(85) 

with the reduced expansion rhn(v) = Y^ilo ^rV 1 and similarly for the other quantities. The 
scaling behavior of the reduced expansions can then be described by 

m R {v) oc (l-/(A)ln(l-^y i/2/3i (l + G(^ln to )) , (86) 
9r(v) oc (l-f(X)\n(l--)) \l + 0(g R \ngR)), (87) 



m R 


= (1-^ 


gn 


= v- 2 g R (v) 


Zr 


= v~ x Z R {v 


Z°r 


= Z°(v), 



Z R {v) ex C(l + 0(g R \ng R )) , (88) 
Z°{v) oc (l-f(\)\n(l--)Y Sl/f3l (l + 0(g R \ng R )), (89) 



which includes improvement. For the purpose of extrapolation and error estimation, auxiliary 
functions similar to eq. (j82"D are being introduced (except in the case of Z R (v), which tends to 
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a constant). From the overlap of the reduced expansion with the auxiliary function one can 
extract the extrapolated value as well as a reliable error estimate. In fact, the proportionality 
factor C = v * rrv R /h& is monotonic and stabilizes for i > 10. With the high-temperature 
expansion up to 13th order, the remainder of the series can be reliably estimated as 

5 = c\h{z)- £> w z4 , (90) 

I i=0 ) z =v/v c 

such that the extrapolated estimate for the (reduced) mass is given by 

fo* R {v) = + (91) 

From the variation in the overlap I obtain an error on 5 and thus on the mass m R . The 
quantities g R and are treated similarly. For Z R , however, I evaluate Z R (v) for L,L — 
1, L — 2, . . . and extract the best estimate by fitting the monotonic behavior of the truncated 
series to a quadratic fall-off in L. 

Results for rriR, g R , Z R , and Z R along the boundary line k = 0.98k c are given in table 6. 
All errors are within 2% and the maximal correlation length is less than 5, compatible 
with the demand that it be less than the average extent of the graphs in the linked cluster 
expansion. As expected, the maximal value for the coupling g R along the boundary line 
occurs in the cr-model limit with g R = 15.8(2), about half the unitarity bound. 

I compare the results from the 13th order calculation with those from the 12th order to 
check for inconsistencies in the error estimates. At 12th order, the critical line is shifted 
to slightly higher k values, but always much less than one standard deviation. The change 
in n c has a negligible effect on m R , g r and Z R , the effect on Z R , however, is noticeable, 
particularly for A < 0.1. This indicates that the error for Z R may be underestimated, but 
this will have no effect on the final result for the triviality bound. 

As a further check, I evaluate the high-temperature series at k — 0.09 in the cx-model limit 
(even though this point is slightly above the boundary line). Here, Monte Carlo data |p~3| 
predict g R = 15.5(1.5) and m R = 0.2075(15). I obtain g R = 15.2(3) and m R = 0.2054(2). 
The Monte Carlo mass is significantly heavier, maybe indicating slight contamination by 
excited states. 

I now turn to the integration of the RG equations fl2"2"|)-([2"5]) in the region 0.98k c < k < n c 
which is done numerically, starting from the results of the high-temperature series at the 
boundary line. The results are displayed in tables 7 — 9 for selected values of A and in the 
range 0.01 < m R < 1.0. The values for k < 0.98n c are from the high-temperature series 
analysis. For the integration I use the RG functions up to three loops including the lattice 
dependent scaling violating terms up to 1-loop order. At this order, the scaling violations 
in the region k > 0.98/t c are always less than 6% of the universal part. Therefore, the 
mass dependent terms are omitted at higher orders. From the tables notice that k —>■ k c 
as m R — > 0, as it should be. Errors are estimated by propagating the errors of the initial 
data and do not include the systematic uncertainties due to the approximation of the RG 
functions. To estimate this systematics, I use the 2-loop approximations instead, and the 
effect for g R in the non-linear limit is demonstrated in fig. 2. Generally, the difference between 
the 2-loop and 3-loop approximations is maximal in the non-linear limit. 
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Table 6: High-temperature expansion results for thr, qr, Zr, and at k — 0.98/t c . 
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3.91597(3) 


0.010910(4) 


0.06 


0.0429388 


0.2787(3) 


1.8492(9) 


3.88140(4) 


0.011060(5) 


0.07 


0.0433295 


0.2776(4) 


2.132(1) 


3.84637(5) 


0.011215(6) 


0.08 


0.0437324 


0.2766(4) 


2.409(2) 


3.81090(7) 


0.011373(7) 


0.09 


0.0441479 


0.2756(4) 


2.680(2) 


3.77500(8) 


0.011535(8) 


0.10 


0.0445763 


0.2747(5) 


2.945(3) 


3.7387(1) 


0.011701(9) 


0.20 


0.0496317 


0.2664(7) 


5.33(1) 


3.3573(3) 


0.01362(2) 


0.30 


0.0561757 


0.2595(8) 


7.34(2) 


2.9653(6) 


0.01608(3) 


0.40 


0.0638780 


0.2535(8) 


9.05(3) 


2.6068(9) 


0.01904(4) 


0.50 


0.0717026 


0.2480(7) 


10.54(4) 


2.321(1) 


0.02221(5) 


0.60 


0.0784481 


0.2429(7) 


11.84(3) 


2.120(1) 


0.02523(6) 


0.70 


0.0834642 


0.2381(6) 


12.98(2) 


1.992(2) 


0.02784(6) 


0.80 


0.0867365 


0.2335(5) 


14.00(1) 


1.915(2) 


0.02999(6) 


0.90 


0.0886387 


0.2291(4) 


14.92(3) 


1.872(2) 


0.03175(5) 


1.00 


0.0897232 


0.2248(2) 


15.8(2) 


1.848(3) 


0.03331(3) 
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Table 7: Results in the symmetric phase at given rriR for A = 0.01. The values for k < 
0.98k c (A) are from the high-temperature expansion, while those for k > 0.98/t c (A) are from 
the integration of the RG equations. 



A rn R 


9r 


Zr 




K 


0.01 1.00 


0.4978807(2) 


4.9641184(1) 


0.00840706(1) 


0.03357427(5) 


0.90 


0.4604509(4) 


4.7752191(2) 


0.00874153(3) 


0.03490241(8) 


0.80 


0.4281371(7) 


4.6061668(3) 


0.00906474(5) 


0.0361834(1) 


0.70 


0.400522(1) 


4.4569599(4) 


0.00937124(9) 


0.0373947(2) 


0.60 


0.377231(2) 


4.3275968(6) 


0.0096553(2) 


0.0385125(3) 


0.50 


0.357933(3) 


4.2180769(8) 


0.0099111(3) 


0.0395125(3) 


0.40 


0.342329(5) 


4.128401(1) 


0.0101333(5) 


0.0403707(4) 


0.30 


0.330145(8) 


4.058578(1) 


0.0103173(8) 


0.0410653(4) 


0.20 


0.321093(9) 


4.008676(1) 


0.0104605(9) 


0.041576(2) 


0.10 


0.314617(9) 


3.978685(1) 


0.0105666(9) 


0.041890(2) 


0.09 


0.314058(8) 


3.976782(1) 


0.0105760(9) 


0.041910(2) 


0.08 


0.313500(8) 


3.975078(1) 


0.0105853(9) 


0.041928(2) 


0.07 


0.312937(8) 


3.973574(1) 


0.0105948(9) 


0.041944(2) 


0.06 


0.312355(8) 


3.972270(1) 


0.0106046(9) 


0.041957(2) 


0.05 


0.311736(8) 


3.971165(1) 


0.0106151(9) 


0.041969(2) 


0.04 


0.311051(8) 


3.970260(1) 


0.0106267(9) 


0.041979(2) 


0.03 


0.310243(8) 


3.969555(1) 


0.0106405(9) 


0.041986(2) 


0.02 


0.309185(8) 


3.969050(1) 


0.0106587(9) 


0.041991(2) 


0.01 


0.307477(8) 


3.968744(1) 


0.0106884(9) 


0.041995(2) 
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Table 8: Same as table 7 but for A = 0.1. 



A 771 z? 

/ \ II UJ-£ 


tf it 


Zn 




K 


0.10 1.00 


4.81569(7) 


4.61008(1) 


0.0091809(1) 


0.0361520(4) 


0.90 


4.4304(1) 


4.43251(2) 


0.0095684(3) 


0.0376002(6) 


0.80 


4.0942(2) 


4.27330(2) 


0.0099494(5) 


0.039001(1) 


0.70 


3.8021(4) 


4.13246(4) 


0.0103196(9) 


0.040330(1) 


0.60 


3.5497(6) 


4.00996(5) 


0.010675(2) 


0.041562(2) 


0.50 


3.332(1) 


3.90583(6) 


0.011011(3) 


0.042669(3) 


0.40 


3.146(2) 


3.82007(8) 


0.011328(5) 


0.043627(3) 


0.30 


2.983(3) 


3.7528(1) 


0.011627(8) 


0.044409(3) 


0.20 


2.833(3) 


3.7043(1) 


0.011926(9) 


0.04499(1) 


0.10 


2.666(3) 


3.67481(9) 


0.01229(1) 


0.04535(1) 


0.09 


2.645(2) 


3.67290(9) 


0.01234(1) 


0.04537(1) 


0.08 


2.623(2) 


3.67117(9) 


0.01239(1) 


0.04539(1) 


0.07 


2.599(2) 


3.66963(9) 


0.01245(1) 


0.04541(1) 


0.06 


2.572(2) 


3.66829(9) 


0.01251(1) 


0.04543(1) 


0.05 


2.541(2) 


3.66712(9) 


0.01259(1) 


0.04544(1) 


0.04 


2.505(2) 


3.66615(9) 


0.01268(1) 


0.04545(1) 


0.03 


2.461(2) 


3.66535(9) 


0.01280(1) 


0.04546(1) 


0.02 


2.402(2) 


3.66472(9) 


0.01296(1) 


0.04547(1) 


0.01 


2.307(2) 


3.66421(9) 


0.01322(1) 


0.04547(1) 
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Table 9: Same as table 7 but for A = 1.0. 



A 771 r> 




Zr 

It 




K 


1.00 1.00 


40.241(3) 


2.3957(2) 


0.0200322(3) 


0.0694724(2) 


0.90 


36.014(6) 


2.2927(3) 


0.0212659(5) 


0.0725736(4) 


0.80 


32.22(1) 


2.1991(5) 


0.022583(1) 


0.0756370(7) 


0.70 


28.81(2) 


2.1149(7) 


0.023996(2) 


0.078616(1) 


0.60 


25.72(3) 


2.040(1) 


0.025525(4) 


0.081454(2) 


0.50 


22.88(6) 


1.975(1) 


0.027200(7) 


0.084092(2) 


0.40 


20.2(1) 


1.920(2) 


0.02908(1) 


0.086459(3) 


0.30 


17.7(2) 


1.875(2) 


0.03129(2) 


0.088481(4) 


0.20 


15.1(3) 


1.841(3) 


0.03408(4) 


0.090073(3) 


0.10 


12.1(1) 


1.819(3) 


0.0383(1) 


0.09113(5) 


0.09 


11.7(1) 


1.818(3) 


0.0389(1) 


0.09121(5) 


0.08 


11.3(1) 


1.816(3) 


0.0395(1) 


0.09127(5) 


0.07 


11.0(1) 


1.815(3) 


0.0402(1) 


0.09133(5) 


0.06 


10.5(1) 


1.814(3) 


0.0411(1) 


0.09139(5) 


0.05 


10.1(1) 


1.812(3) 


0.0420(1) 


0.09143(5) 


0.04 


9.59(9) 


1.811(3) 


0.0432(2) 


0.09147(5) 


0.03 


9.02(8) 


1.810(3) 


0.0446(2) 


0.09150(5) 


0.02 


8.31(7) 


1.809(2) 


0.0466(2) 


0.09152(5) 


0.01 


7.33(5) 


1.808(2) 


0.0498(3) 


0.09154(5) 
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Figure 2: Comparison of the solutions of the RG equations at two and three loops for gR at 
A = 1. The dashed line is the result from the high-temperature expansion, while the solid 
line is the continuation obtained by integrating the 3-loop approximation of the /3-funxction. 
The dotted curve corresponds to the 2-loop approximation. It is evident that the deviations 
in the symmetric phase are of the same size as the series expansion error in gR calculated at 
the boundary line. 



4.2 Continuation of the analysis to the broken symmetry phase 

As indicated in section 2, to continue the solution of the theory into the broken phase one 
needs to evaluate the integration constants Ci, C2, and C3 at the transition. These constants 
are defined through the scaling behavior of ibr, Zr, and Z® in eqs. (f|), ( ^6l) and (|27| ) and can 
be computed by numerically integrating the RG equations down to small mjj, for instance 
tur = 10~ 8 , using the 3-loop expressions. In this region InCj is linear in gR and the limit 
gn — > can be extracted. The results for some values of A are presented in table 10. 

Given the integration constants on the high-temperature side of the transition, the con- 
stants C[ in the broken phase are fixed through eq. PSf ). tur, Zr, and Z^ can then be 
calculated using the broken phase scaling laws. I choose Qr = 10~ 8 — small enough for 
order g\ corrections to be negligible — and continue integrating the RG equations up to 
itlr = 0.5, again using the 3-loop expressions for the RG functions including scaling violating 
terms up to 1-loop. The results at rriR = 0.5 are displayed in table 11 for selected values of 
A. Results along the lines of constant A are shown in tables 12 — 14. Errors are estimated 
by propagating the uncertainties in the integration constants only. The renormalized cou- 
pling is maximal in the non-linear limit where it reaches gR = 18(1) at tur = 0.5, about 
2/3 the unitarity bound. The assumption of applicability of perturbative RG is justified by 
comparing the 3-loop with the 2-loop results: the changes in the region rriR < 0.5 are not 
unreasonable. 



In fig. 3 I compare the results for tur/ ' with recent Monte Carlo data [111]. Both methods 
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Table 10: Values for InC; at the critical line. 



A 


InCi 




lnC 2 


lnC 3 


0.01 


126.429(4) 


1 


3782607(3) 


-5.1270(1) 


0.02 


64.356(5) 


1 


370040(1) 


-4.7730(2) 


0.03 


43.563(7) 


1 


361631(2) 


-4.5625(3) 


0.04 


33.116(8) 


1 


353031(4) 


-4.4106(5) 


0.05 


26.82(1) 


1 


344237(6) 


-4.2907(6) 


0.06 


22.60(1) 


1 


335249(9) 


-4.1909(7) 


0.07 


19.57(1) 


1 


32606(1) 


-4.1049(9) 


0.08 


17.29(1) 


1 


31668(1) 


-4.029(1) 


0.09 


15.51(1) 


1 


30710(2) 


-3.961(1) 


0.10 


14.07(2) 


1 


29731(2) 


-3.898(1) 


0.20 


7.50(2) 


1 


18865(8) 


-3.441(3) 


0.30 


5.22(2) 


1 


0636(2) 


-3.110(4) 


0.40 


4.04(2) 





9340(3) 


-2.831(4) 


0.50 


3.32(1) 





8173(4) 


-2.596(4) 


0.60 


2.82(1) 





7263(6) 


-2.407(4) 


0.70 


2.460(7) 





6633(8) 


-2.259(3) 


0.80 


2.184(4) 





624(1) 


-2.144(3) 


0.90 


1.963(7) 





601(1) 


-2.052(3) 


1.00 


1.78(4) 





588(1) 


-1.975(9) 
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Figure 3: Ratio tur/ fir versus itir in lattice units. The solid line is from integrating the 
RG equations. The dashed line gives the propagated error and the dotted line is the 2-loop 



result. The data points are from a Monte Carlo simulation [13j 



are in fair agreement for uir < 0.5. For larger masses, however, one leaves the perturbative 
regime and the integration of the perturbative RG functions cannot be trusted anymore. 
Also, the L — > oo extrapolation of the Monte Carlo data may be spoiled by cx-particle decay, 
perhaps leading to the systematic difference evident in fig. 3. At A/itlr = 21 obtain for an 
upper bound on the Higgs mass 

m R /U < 2-46 ± 0.02 HTE ± 0.08 PT , (92) 

where the first error is from the extrapolation of the high-temperature expansion and the 
second error is from the perturbative expansion of the RG functions. The latter is estimated 
by comparing the 2-loop and the 3-loop results. The above numbers translate to about 
600 ± 5hte ± 20pt GeV. The Monte Carlo result is 2.3(2) Hl3| , somewhat below the number 
from the high-temperature expansion. The agreement is better in the second reference of [|T3" 



where in a similar plot the RG integration is based on Monte Carlo data in the symmetric 
phase. This indicates that the systematic difference in fig. 3 is mainly due to the expansion 
itself and not to the perturbative RG integration. 

The F4 results compare to vtlrI < 2.6(3) (or 640(70) GeV) from the high-temperature 
expansion on a hypercubic lattice || 3. Here Monte Carlo simulations give 2.7(1) [ |T0"|| . 
slightly above the semi-analytical result. Interestingly, the hypercubic results are in reverse 
order. 

Even though the maximal cutoff on the _F 4 lattice seems to be about 6% lower than 
LW's hypercubic result [|J , this statement by itself is not meaningful. One must not directly 
compare the two cutoff's as they have different meanings on different lattices. This is evident 
in the lattice dependence of the scaling violations. Rather, one should consider the cutoff 
dependence of physical quantities, such as the pion — pion 90 degree scattering cross section. 
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Table 11: Results in the broken phase at m R = 0.5 from integrating the RG equations. 



A 


K 


9r 


Z R 


7° 
^R 


0.01 


0.043384(1) 


0.290944(8) 


3.841565(1) 


0.010567(1) 


0.02 


0.043739(3) 


0.57660(5) 


3.810283(4) 


0.010709(3) 


0.03 


0.044106(4) 


0.8573(1) 


3.778535(9) 


0.010855(4) 


0.04 


0.044484(5) 


1.1335(3) 


3.74633(2) 


0.011004(6) 


0.05 


0.044874(6) 


1.4053(5) 


3.71368(2) 


0.011156(8) 


0.06 


0.045276(8) 


1.6731(8) 


3.68059(3) 


0.01131(1) 


0.07 


0.045691(9) 


1.937(1) 


3.64707(4) 


0.01147(1) 


0.08 


0.04612(1) 


2.197(2) 


3.61313(6) 


0.01164(2) 


0.09 


0.04656(1) 


2.454(2) 


3.57879(7) 


0.01180(2) 


0.10 


0.04701(1) 


2.708(3) 


3.54407(8) 


0.01198(2) 


0.20 


0.05238(2) 


5.08(1) 


3.1800(3) 


0.01397(5) 


0.30 


0.05933(2) 


7.21(3) 


2.8068(5) 


0.01652(9) 


0.40 


0.06752(3) 


9.15(4) 


2.4659(8) 


0.0196(1) 


0.50 


0.07585(3) 


10.92(4) 


2.195(1) 


0.0229(1) 


0.60 


0.08305(4) 


12.55(4) 


2.004(1) 


0.0261(1) 


0.70 


0.08844(4) 


14.05(3) 


1.882(2) 


0.0288(1) 


0.80 


0.09199(5) 


15.46(2) 


1.809(2) 


0.0310(1) 


0.90 


0.09409(4) 


16.79(4) 


1.768(2) 


0.0329(1) 


1.00 


0.09534(2) 


18.1(3) 


1.745(2) 


0.0345(6) 



Table 12: Results in the broken phase at given m R for A = 0.01. 



A m R 


9r 


Zr 




K 


0.01 0.01 


0.306186(8) 


3.967302(1) 


0.010722(1) 


0.041998(2) 


0.02 


0.307806(8) 


3.967149(1) 


0.010693(1) 


0.042000(2) 


0.03 


0.308731(8) 


3.966898(1) 


0.010676(1) 


0.042002(2) 


0.04 


0.309357(8) 


3.966549(1) 


0.010664(1) 


0.042006(2) 


0.05 


0.309810(8) 


3.966101(1) 


0.010655(1) 


0.042011(2) 


0.06 


0.310147(8) 


3.965555(1) 


0.010647(1) 


0.042017(2) 


0.07 


0.310396(8) 


3.964909(1) 


0.010641(1) 


0.042023(2) 


0.08 


0.310576(8) 


3.964165(1) 


0.010635(1) 


0.042031(2) 


0.09 


0.310698(8) 


3.963322(1) 


0.010630(1) 


0.042040(2) 


0.10 


0.310770(8) 


3.962380(1) 


0.010626(1) 


0.042051(2) 


0.20 


0.309520(8) 


3.947501(1) 


0.010597(1) 


0.042211(2) 


0.30 


0.305563(8) 


3.922586(1) 


0.010580(1) 


0.042481(2) 


0.40 


0.299318(8) 


3.887393(1) 


0.010570(1) 


0.042869(2) 


0.50 


0.290944(8) 


3.841565(1) 


0.010567(1) 


0.043384(1) 
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Table 13: Same as table 12 but for A = 0.1. 





!JR 


n R 


^R 




10 01 


9 917(9) 




01 3^4(9) 


U .\J L ±0 L ±cf \ ± J 


02 


2 32^(2) 


3 65471 (9) 


01 328(2) 


D D4549(1 ) 


03 






01 31 9(9) 


fl 04^40(1 1 


04 

U.Ut: 


9 A90(9) 




01 301 (9) 


04^ ^(1(1 ) 


0.05 


2 453(2) 


3 65349(9) 


o 01293(2) 


04550(1) 


0.06 


2.479(2) 


3.65293(9) 


0.01286(2) 


0.04551(1) 


0.07 


2.502(2) 


3.65230(9) 


0.01279(2) 


0.04552(1) 


0.08 


2.522(2) 


3.65158(9) 


0.01274(2) 


0.04553(1) 


0.09 


2.540(2) 


3.65078(9) 


0.01269(2) 


0.04554(1) 


0.10 


2.556(2) 


3.64989(9) 


0.01265(2) 


0.04555(1) 


0.20 


2.658(3) 


3.63640(9) 


0.01237(2) 


0.04573(1) 


0.30 


2.705(3) 


3.61440(9) 


0.01220(2) 


0.04603(1) 


0.40 


2.719(3) 


3.58373(9) 


0.01207(2) 


0.04646(1) 


0.50 


2.708(3) 


3.54407(8) 


0.01198(2) 


0.04701(1) 



Table 14: Same as table 12 but for A = 1.0. 



A 1TLR 


9r 


Zr 




K 


1.00 0.01 


6.64(4) 


1.793(2) 


0.0538(7) 


0.09156(5) 


0.02 


7.43(5) 


1.792(2) 


0.0509(6) 


0.09156(5) 


0.03 


7.99(6) 


1.792(2) 


0.0491(6) 


0.09158(5) 


0.04 


8.43(7) 


1.791(2) 


0.0478(6) 


0.09159(5) 


0.05 


8.81(8) 


1.790(2) 


0.0468(6) 


0.09161(5) 


0.06 


9.14(8) 


1.790(2) 


0.0460(6) 


0.09163(5) 


0.07 


9.45(9) 


1.789(2) 


0.0452(6) 


0.09166(5) 


0.08 


9.73(9) 


1.789(2) 


0.0446(6) 


0.09169(5) 


0.09 


10.0(1) 


1.788(2) 


0.0440(6) 


0.09172(5) 


0.10 


10.2(1) 


1.787(2) 


0.0435(6) 


0.09176(5) 


0.20 


12.3(1) 


1.780(2) 


0.0400(6) 


0.09230(4) 


0.30 


14.1(2) 


1.770(2) 


0.0377(6) 


0.09312(2) 


0.40 


16.0(2) 


1.758(2) 


0.0360(6) 


0.09415(2) 


0.50 


18.1(3) 


1.745(2) 


0.0345(6) 


0.09534(2) 
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By demanding that the deviation 



(da/d6) laW " 
(da/d&f 



i r i i n\ continuum 

(93) 



be less than (say) 0.3% for all energies up to some center-of-mass energy, for instance up 
to 2m#, one defines a physically relevant bound. On the lattice, 5 < 0.3% translates to 
mji pa 0.5, whereas on the hypercubic lattice the same requirement gives mjj px 0.22 [12|. 



From ||, this value corresponds to m R / f n = 2.2 or roughly 540 GeV, which is actually about 
10% below the F 4 bound. 

To summarize, in the broken phase I find that for A/m^ > 2 the coupling never exceeds 
2/3 of the triviality bound (as in the hypercubic case, see ||), indicating that there is in 
fact not much room for a heavy Higgs — if one insists on keeping the naive lattice action. 
A Monte Carlo study including dimension six operators parameterizing the leading order 
(A -2 ) cutoff effects has shown that the bound may be as high as 710(60) GeV [I4"|. A high- 



temperature expansion for the improved action case would be desirable but poses quite a 
complicated task. 

Going from 12th to 13th order in the expansion quadruples the number of contributing 
graphs and one reaches the limits of todays computational resources. Furthermore, at this 
order, the uncertainty due to the perturbative RG functions is by far larger than the uncer- 
tainty due to the series expansion. Hence it does not pay to attempt a 14th order calculation 
without simultaneous improvement in the perturbative results for the RG functions. 

Finally, I would like to remind the reader that the solution of the 4 theory as presented 
here rests on the assumption that perturbation theory can be applied for sufficiently small 
couplings (thereby implicitly assuming triviality), and relies on the fact that the perturbative 
region overlaps with the region in the symmetric phase where the linked cluster expansion 
can be evaluated. These issues are discussed in detail in LW's papers. 

The advantages of the semi-analytical approach over Monte Carlo methods are evident: 
the lack of finite size effects and the absence of statistical errors. Furthermore, one does not 
have to worry about the cr-particle decay which complicates the extraction of m a on larger 
lattices in Monte Carlo simulations. However, with respect to fig. 3 it is not clear from which 
source the differences arise. All one can say is that the comparison gives a good estimate of 
systematic errors which are very different in both methods and presumably do not alter the 
results in the same direction. 



Acknowledgements 

I would like to thank my former advisor, Urs M. Heller for invaluable support. Thanks 
also go to Pavlos Vranas and Herbert Neuberger. This project was supported in part by 
the Department of Energy under contract number DE-FC05-85 ER250000 and the Super- 
computer Computations Research Institute and The Florida State University Department 
of Physics through use of their facilities. The linked cluster expansion was performed on 
the cluster of IBM RS6000 workstations at SCRI. The final stage of this project was car- 
ried out while the author was at Columbia University, supported under contract number 
DE-FG02-92 ER40699. 



31 



Appendix 

This appendix gives the (slightly corrected) coefficients for the RG function (cf. |5|, |12|). 
Writing 

oo 

P(0,9r) = 9rY,Pv9r, (94) 

v=l 

oo 

7(0,<ta) = £7„<&, (95) 

u=l 
oo 

6(0, g R ) = J2 6 »9r, (96) 



the universal coefficients are given by 



Pi = |(iV + 8)(167r 2 )- 1 ) (97) 

ft = -^(3iV+14)(167r 2 )- 2 , (98) 
/3 3 = -^-|61iV 2 + 1782A^ + 5744 + 192C(3)(5iV + 22) 

-2A;(iV + 8)(13A^ + 62))(167r 2 )- 3 , (99) 



71 = 0, (100) 

72 = ±:(N + 2)(l(rf)-\ ( 101 ) 
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73 = -^(iV + 2)(iV + 8)(l-2A;)(167r 2 )- 3 , (102) 

51 = -i(iV + 2)(167r 2 )- 1 , (103) 

5 2 = l(iV + 2)(167r 2 r 2 , (104) 

lo 



6 



?> 



2 ^(iV + 2)|31iV + 221 + 6k(N + l)|(16vr 2 )- 3 , (105) 



where 

k = 1.62520965, (106) 
C(3) = 1.20205690, (107) 

the latter being the Riemann zeta function. The non-universal scaling violating terms up to 
one loop are 



P(m R ,g R ) = gR^2u y {m R )g R , (108) 

oo 

l(m R ,g R ) = ^v u {m R )g R} (109) 

oo 

6{m R ,g R ) = J2 w "( m R)9R, (HO) 
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with the coefficients 



4m| 



4 + mfj 

«i = jU {{N + 8)(4 + m^)J 3 (m^) - QJ 2 {m 2 R ) 

+ Mml)} , (112) 

1 . , 

= 7% , (113) 



Vl 



= + 2) Mo { J 2 (m|) - —Up JiK)} , (114) 

1 / X 

w o = -7«o , (115) 



Wl = -2v 1 -^(N + 2)m 2 R {(A + m 2 R )J 3 (m 2 R )-2J 2 (m 2 R ) 



The lattice integrals J p (m R ) are defined as 



py " lR) = L [g(k) + mlV (U7) 



JJm 2 



and are evaluated numerically. 

In the broken phase, the universal coefficients fli, /3 2 , 7i, and <5i take the same form as 
in the symmetric phase; the differing coefficients are given by 

(3 3 = (-— N 3 - 0.6038002iV 2 + 6.006124iV+ 11.1064ll(167r 2 )- 3 , (118) 
I 108 J 

72 = -^(16tt 2 )- 2 , (119) 



73 = j-0.0452317896iV + 0.211712032 j(167r 2 )- 3 , (120) 

5 2 = ^j(5 + 2v / 37r)iV + 16-2v / 37rj(167r 2 )- 2 , (121) 

5 3 = (^-iV 3 + 1.2932334Ar 2 -4.320590Ar + 1.620019l(167r 2 )- 3 . (122) 

1 108 J 

For the function e, the power series is written as 

<0,9r) = f><&, (123) 

with the universal coefficients 

eo = 1, (124) 
e = i|-iV + 3-2v / 37r|(167r 2 )- 1 , (125) 

e = |o.503662A^ + 2.487506|(167r 2 )^ 2 . (126) 
33 



For completeness, I quote the scaling violating terms (taken from El) up to terms of order 



m R In m R 



m R smhm R m R smhm R 

u o = 2 - ■ , 2l o r TJ~i ' ( l27 > 

2 sinh \m R 2 — sinh \m R 



2 f , 2 »r / 7 ^ 47r2 167r2 
Mi = ft + m R ^-2\nm R + N I — + —r - — n — r 2 

13 7T 2vr 2 112vr 2 . 704tt 2 1 2 w 



+ 36 - 7! + ^ ^ n + ^ r2 | (167F 

+ %fllnm^ (128) 
1 hir sinh m R 
4 2 — sinh 2 



IIOROlllLl IIV R / 10 M 

= -7 ^ — 2 ; ~ , ( 129 ) 



2 fl /lT N1 2 /vr 2 4tt 2 \ 2tt 2 8tt 2 

Wl = m H^( iV + 2 ) lnm ^ + iV (y ro + — ri ) +ir ro + ~ n 

64vr 2 3^ - 4 
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r 2 + ^-^{(lGvr 2 )- 1 + 0(m R \nm R ) , (130) 

m R sinh m R sinh 2 £m B 

= r-T4i > ( 131 ) 

4 — sinh im/j 

2 fl ar nNl 2 Ar /vr 2 16tt 2 16tt 2 7\ 
Wl = 5 1 + m|| — (iV + 2) In m R + Nl—r + -_ ri + —r 2 - — J 

2tt 2 32vr 2 224vr 2 5 V^rl,.,,, 2 n-i 

" r 7^ r l n r 2 - 77 + "77" H167T ) 



3 u 9 9 18 36 

+ 0(m R \nm 2 R ), (132) 

where the constants tq, r±, r 2 are given by 

r = 0.13823047, (133) 
ri = -0.04029906 , (134) 
r 2 = 0.00763417. (135) 

For e, the m^-dependent expansion is written as 



£(m R ,g R ) = ^x u (m R )g R , (136) 



u=0 



with the coefficients 



sinh 77^ 2 + sinh 2 %m R 
^0 = 77727 ' ( 137 ) 



m R 2 - sinh im R 



2 f, 2 »,/ 1 27T 2 

= ei + m^jlnm^ + N I -— - yr + — n + — r 2 



8tt 



2 



7T 2 207T 2 1127T 2 11 7^. 

~y r ° + — ri " — r2 -24 + ^6- Kl67r 



2\ 



+0(m^lnm 2 J ). (138) 
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